ASYMPTOTICS FOR HESSENBERG MATRICES FOR THE 
BERGMAN SHIFT OPERATOR ON JORDAN REGIONS 



EDWARD B. SAFF AND NIKOS STYLIANOPOULOS 

Abstract. Let G be a bounded Jordan domain in the complex plane. 
The Bergman polynomials {p n }^L of G are the orthonormal polynomi- 
als with respect to the area measure over G. They are uniquely defined 
by the entries of an infinite upper Hessenberg matrix M. This ma- 
trix represents the Bergman shift operator of G. The main purpose of 
the paper is to describe and analyze a close relation between M and 
the Toeplitz matrix with symbol the normalized conformal map of the 
exterior of the unit circle onto the complement of G. Our results are 
based on the strong asymptotics of p n . As an application, we describe 
and analyze an algorithm for recovering the shape of G from its area 
moments. 



1. Introduction 

Let G be a bounded simply-connected domain in the complex plane C, 
whose boundary T := dG is a Jordan curve and let {p n }%Lo denote the 
sequence of Bergman polynomials of G. This is the unique sequence of 
polynomials 

p n (z) = X n z n + ■ ■ ■ , A n >0, 71 = 0,1,2,..., (1.1) 
that are orthonormal with respect to the inner product 



</,<?> == / f(z)g(z)dA(z), 

JG 

where dA stands for the area measure. We denote by L 2 a {G) the Hilbert 
space of all functions / analytic in G for which 

WfWmo ■= (/,/) 1/2 <oo, 

and recall (cf. [6]) that the polynomials {p n }'^Lo form a complete orthonor- 
mal system for L 2 a {G). 
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Let Q, := C\G denote the complement of G in C and let <I> denote the 
conformal map Q — > A := {w : \w\ > 1}, normalized so that near infinity 



Finally, let ^ 



with 



= : A 



iz + 70 H 1 — H , 7 > 0. 

z Z L 



(1.2) 



6u> 



O denote the inverse conformal map. Then 
62 



&o + — ■+ 



+ 



M > 1, 



w 



iv- 



(1.3) 
(1.4) 



& = I/7 = cap(r), 
where cap(r) denotes the (logarithmic) capacity of T. 

On L^(G) we consider the multiplication by z operator (also known as 
the Bergman shift operator) A4 : / — > zf. Note that A4 defines a bounded, 
noncompact, linear operator on L^(G) and that 

a ess (M) = T; (1.5) 

see [1], where we use a ess (L) to denote the essential spectrum of a bounded 
linear operator L; that is, the set of all A £ C for which L — XI is not a 
Fredholm operator. For the operators we consider, the essential spectrum is 
the same as the continuous spectrum. 

We also consider the matrix representation of Ai in terms of the orthonor- 
mal basis {p n }^=o- This induces the upper Hessenberg matrix 



M 



'on 


&01 


bo2 


&03 


^04 


»10 


&11 


b\2 


&13 


614 





021 


b22 


^23 


&24 








&32 


^33 


&34 











&43 


644 



(1.6) 



where 



h,j = (zpj,Pk), k>0, j > 0. 



Note that b^ j = for k> j + 2 and that 



n+l 



Zp n (z) = y^&fc,nPfcQ 



fc=0 



In particular, 



It follows 



An 



+ 1 



>0, n = 0,1,.... 



b n+1<n p n+ i(z) = Zp n {z) - y^6fc,nPfc(^ 



(1.7) 
(1.8) 

(1.9) 
(1.10) 



k=0 



and, hence, the entries of M define uniquely the sequence of Bergman poly- 
nomials of G, in the sense that p n +i, n = 0, 1, . . ., can be computed recur- 
sively from (1.10). 
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It is shown in [11] and [9] (see also [12, Thm 2.4]) that except for some 
trivial cases, the matrix (1.6) is not banded; i.e., the p n 's do not satisfy 
a recurrence relation of bounded length. It is also well-known that the 
eigenvalues of the n x n principal submatrix of M coincide with the zeros 

Of p n {z). 

Our goal is to investigate the asymptotic behavior of the entries in the 
matrix M. In particular, we show that if the boundary of G is piecewise 
analytic without cusps, then all the diagonals (sub, super and main) have 
limits which are the coefficients of the Laurent expansion (1.3) of the inverse 
conformal map \P: 

lim b n+ln = b and lim b n _ kn = b k , k = 0,l, (1-11) 

n— >oo n— yoo 

A potential application of (1.11) is in the area of geometric tomography, 
where the following inverse problem arises: Given a finite number of complex 
moments 

fi kj :=(z k ,zJ)= [ z k zidA(z), k, j = 0, 1, . . . , (1.12) 
JG 

how can one approximate the region G that generated these moments? Re- 
garding existence and uniqueness, we note a result of Davis and Pollak [4] 
stating that the infinite matrix \p m ,k\m k=o defines uniquely the curve T. By 
utilizing the given moments to compute Bergman polynomials, and thereby 
a principal submatrix of M, the subdiagonals of the submatrix will provide 
an approximation to the Laurent coefficients of the mapping of the unit 
circumference onto the boundary of G. We will discuss this procedure in 
Section 3. 

We note that there is a one-to-one correspondence between the complex 
moments (1.12) and the real moments 



T mn ■= / x m y n dxdy, m,n = 0,1, (1-13) 

JG 

Namely, 

J2J2i m - j i n ~ k ( m r \r j+kim+n _ j _ k , i :=«/=!, (1.14) 

„■— n i — n \ J / \ J 



m n 

l-l"m,n — 

j=0 k=0 

or, in the inverse direction, 



Tm,n = (^EEffl (1-15) 

j= ok=o VJ/ V fc / 

see [4]. Thus, the moments in (1.12) will uniquely determine the moments 
in (1.13) and vice-versa. 

The Faber polynomials {-F n }^L of G are defined as the polynomial part 
of the expansion of ^"(z), n = 0, 1, . . ., near infinity, that is, 

<Z> n (z) = F n (z)-E n (z), (1.16) 
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F n (z) = 7 n z n + 



and E n (z) = O 



z — > oo. 



(1.17) 



The Faber polynomial of the 2nd kind, G n (z), is defined as the polynomial 
part of $ n (z)$'(z), that is, 

G n (z) = * n (z)&{z) - H n (z), z € n, 

where 

1 



G n (z)= 7 n+l z n + 



and H„ (z) 



O 



z — > oo. 



(1.18) 
(1.19) 



It follows immediately from (1.16) and (1.18) that 

G n (z) = %^ and H n (z) = ^±^-. (1.20) 
n + 1 71 + 1 

It is well-known that the Faber polynomials of the 2nd kind satisfy the 
following recurrence relation (see [5, p. 52]): 

n 

zG n (z) = bG n+1 (z) + Y / b j G n - j (z), G (z) = b. (1.21) 

3=0 

Consider now the Toeplitz (and upper Hessenberg) matrix defined by 
the continuous function on T := {w : \w\ = 1}, that is, 



T* := 



b 




h 


h 


b 4 


b 


b 


bi 


h 


h 





b 


h 


bi 


b 2 








b 


bo 


h 











b 


bo 



(1.22) 



It follows from (1.21) that the eigenvalues of the nxn principal submatrix of 
T$ coincide with the zeros of G n (z); see also [15]. This is a relation similar 
to the one connecting the upper Hessenberg matrix M with the Bergman 
polynomials {p n }%L - 

In [13, §7.8] it is shown that if T is piecewise analytic without cusps, then 



< 



ci(T)- 



n E N, 



(1.23) 



where urn (0 < oj < 2) is the smallest exterior angle of T. (Hereafter, we 
use Cfc(r), k = 1,2,..., to denote a non-negative constant that depends 
only on T.) Therefore, in this case, the symbol ^ of the Toeplitz matrix Ts^ 
belongs to the Wiener algebra, which leads to the conclusion that T<j, defines 
a bounded linear operator on the Hilbert space I 2 and that 

a ess (T*) = F; (1.24) 

see e.g. [2, p. 1-10]. 
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We end this section by noting a result, regarding a property of H n , that 
we are going to use in Section 5. A proof can be found in [13, Lem. 2.1]. 

Lemma 1.1. For any n G N, H n is analytic and square integrable in Vt. 

2. Main results 

In this section we state and discuss our main results. Their proofs are 
given in Section 5. Section 3 contains applications of our results to the 
recovery of planar regions. 

From (1.5) and (1.24) it follows that 

a ess (M) = a ess (T^). (2.1) 

The next theorem shows that the connection between the matrices M and 
is much more substantial. 

Theorem 2.1. Assume that T is piecewise analytic without cusps. Then, it 
holds as n — )• oo, 

2 b n+hn = b + o(-), (2.2) 
V lit J 

and for k > 0, 

V n+l ' n ' - Vv^, 

where O depends on k but not on n. (See (5.17) for more precise estimates.) 

Improvements in the order of convergence occur in cases when V is smooth. 
In order to state the corresponding results we need to introduce the smooth- 
ness class C(q, a) of Jordan curves. We say that T belongs to C(q,a), q E N, 
if r is defined by z = g(s), where s denotes arclength, with <jrW g Lip a, 
for some < a < 1. Then both $ and VP := are q times continuously 
differentiable in O \ {oo} and A \ {oo} respectively, with (fr( Q ^ and in 
Lip a: see, e.g., [14, p. 5]. 

Theorem 2.2. Assume that T € C(p + l,a), with p + a > 1/2. Then, it 
holds as n — > oo, 



n+l \ n 



n — k + 1 , , ^ / 1 

b n -k,n = b k + O [-= , (2.3) 



1 &»-*,» = 4 * + ° l^rr], (2.5) 



and for k > 0, 



n+l ' V reP+ °, 

where O depends on k but not on n. (See (5.29) for more precise estimates.) 

For the case of an analytic boundary V further improved asymptotic re- 
sults can be obtained. To state these results we need to introduce some 
notation. For an analytic curve Y the mapping \& can be analytically con- 
tinued as a conformal map to the exterior of some disk {w : \w\ < g}, where 
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< q < 1. We denote by L a the image of the circle {w : \w\ = a} under the 
map ty. In other words, 

L a := {z € C : |$(z)| = a}. 

Theorem 2.3. * Assume that the boundary T is analytic and let g < 1 be 
the smallest index for which <3? is conformal in the exterior of L g . Then, it 
holds as n — > oo, 

n + 1 

and /or k > 0, 



b n+1 , n = b + 0( e 2n ), (2.6) 



n-k + 1 



b n -k,n = b k + 0{y/n\ogng n ), (2.7) 

n + 1 

where O depends on k but not on n. (See (5.39) for more precise estimates.) 

In the converse direction we have: 

Theorem 2.4. Assume that V is a Jordan curve without zero interior an- 
gles. If 

l/n 



lim sup 

n— >oo 



n + 2 h h 

— rO n +l ,n — 

n+1 



< 1, (2. 



then r is analytic. 



The following example shows that the inverse statement does not make 
sense for the main diagonal of M. 

Example 2.1. Consider the case where the domain G has m-fold rotational 
symmetry about the origin, for some m > 2. 

This means that Q l27T / m z £ O, whenever z £ £1. Then, it is easy to see 
that 

&0 = and b n ^ n = 0, n > m. (2-9) 
Indeed, by using symmetry arguments it follows 

^(e i2w/m w) = e i27r/m y(w), wen, (2.10) 

and for n = km + j, with j = 0, 1, . . . , m — 1, 

Pn (z) = zUk(z m ), deg(g fc ) = k. (2.11) 

The first relation in (2.9) follows at once from (2.10). For the second relation 
in (2.9), observe that (2.11) implies for n > m that 

p n {z) = \ n z n + 0{z n - m ), 

which, in turn, yields (z n+ ,p n ) = and therefore b n ^ n = (zp n ,p n ) = 0. 



*This theorem, along with a sketch of its proof given in Section 5.3, was presented by 
the first author at the Joint Meeting of the AMS and MAA in Phoenix, January 2004. 
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3. A RECOVERY ALGORITHM 
Reconstruction Algorithm 

1. Start with a finite set of complex moments fikj, k,j = 0,1,..., n; see 
(1.12), or, equivalently from a finite set of real moments t^j, k,j = 
0,1,..., n; see (1.13). 

2. Use the Arnoldi version of the Gram-Schmidt ( GS) process, in the way 
indicated in [13, §7.4] , to construct the Bergman polynomials {pk}k=o 
from the moments fj,f.j, k,j = 0,1,..., n. This involves at the k-step 
the orthonormalization of the set {po,pi, ■ ■ ■ ,Pk-l, zp^-i}, rather than 
the set of monomials {1, z, . . . , z k ~ 1 , z k }, as in the conventional GS. This 
process, in particular, yields the inner products 

h,j = (zpj,Pk), 3=0, l,...,n, k = 0,...,j + 1. 

3. Choose a number m, 1 < m < n, and set 

i \ In + 2 r n ) In — k + 1 , 

o v ; := \ — —b n+ i. n , bi':=\ — — b n - k . n , k = 0, 1, ...,m. 

Vn+1 V^+l 

(3.1) 

(See Theorem 3.1 and Remark 3.2 below, for a suitable choice ofm.) 

4. Form 

^ ) H:=^ + ^ ) + V + --- + ^- (3 ' 2) 

5. Approximate V by Tm , where 

r&> :=¥&)(«;), a)£l (3.3) 

Remark 3.1. We refer to [13, §7.4] for a discussion regarding the stability 
properties of the Arnoldi GS. In particular, we note that the Arnoldi GS does 
not suffer from the severe ill- conditioning associated with the conventional 
GS as reported, for instance, by theoretical and numerical evidence in [10]. 

The following result justifies the use of the algorithm for analytic curves. 

Theorem 3.1. Assume that T is analytic, and let g < 1 be the smallest 
index for which $ is conformal in the exterior of L Q . Set n = 2m. Then, 
for any \w\ > 1 it holds that 

|*H - V$(w)\ < ci(rVmlogmf? m + c 2 {T)\w\ Q im , (3.4) 
where the constants c\(T) and C2(T) depend on T only. 

Remark 3.2. Similar estimates, as in the above theorem, can be obtained 
for the case where T is piecewise analytic without cusps. However, these 
estimates are too pessimistic compared with actual numerical evidence; see 
Figure 2 below. We were only able to rigorously show that for an uniform 
error of order 0(1/ \/m) we require the computation of the orthonormal poly- 
nomials up to degree m 4+a; , where uj-k is the smallest exterior angle ofT. 
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Figure 1. Recovery of an ellipse, with n = 10 (left) and 
n = 20 (right). 



For applications to the 2D image reconstruction arising from tomographic 
data we refer to [8]. Here we highlight the performance of the reconstruc- 
tion algorithm by applying it to the recovery of three curves, coming from 
different classes: an analytic curve, one curve with corners and one curve 
with cusps. For providing matter for comparison with the reconstruction al- 
gorithm of [13, §7.7] we have chosen to present results for exactly the same 
curves as in [13]. We note that the reconstruction algorithm of [13] is based 
on approximating first the exterior conformal mapping w = &(z) in terms of 
the ratio p n+ \{z) / p n {z) , cf. the estimates (5.3)-(5.4) below, and then on in- 
verting the so-formed Laurent series in order to compute an approximation 
of the inverse map z = ^f(w). 

In each case we start by computing a finite set of complex moments (1.12) 
up to degree n, and then follow the steps 2—5 of the algorithm, taking 
m = n/2. In all three examples the complex moments are known explicitly. 
All computations were carried out on a desktop PC, using the computing 
environment MAPLE. 

(n) 

In Figures 1-3 we depict the computed approximation Tm against the 
original curve T. The presented plots indicate that the above reconstruction 
algorithm constitutes a valid method for recovering a shape from its partial 
moments. Even in the cusped case, pictured in Figure 3, the fitting is 
remarkably close, despite the low degree of the moment matrix used. 

In Figure 1 we illustrate the reconstruction of an ellipse, where, with the 
notation of Theorem 3.1, g = 1/3. 

In Figure 2 we reconstruct a square by using the complex moments up 
to the degree 16. We have chosen n = 16, so that the result can be com- 
pared with the recovery of a square, as shown on page 1067 of [7], obtained 
using the Exponential Transform Algorithm. This is another reconstruction 
algorithm based on moments. 

In order to show that the proposed reconstruction algorithm works equally 
well for domains where the results of neither Theorem 3.1 nor that of Re- 
mark 3.2 apply, we use it for the recovery of the boundary of the 3-cusped 
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Figure 2. Recovery of a square, with n = 16. 




Figure 3. Recovery of a 3-cusped hypocycloid, with n = 20 
(left) and n = 30 (right). 

hypocycloid defined by 

r := {z = V(w) = w + weT}. 

The application of the algorithm with n = 20 and n = 30 is depicted in 
Figure 3. 

Comparing the performance of the above algorithm with that of [13] for 
the cases of the ellipse and the hypocycloid, it appears that the latter algo- 
rithm performs slightly better. On the other hand, both algorithms perform 
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better than the reconstruction algorithms of [7] for the case of the square. 
More definitive comparisons will require further experimentation and anal- 
ysis of all three reconstruction algorithms. 

4. Numerical Results 

In this section we employ the first three steps of the reconstruction al- 
gorithm in order to present numerical results that illustrate the order of 
convergence in (2.2) and (2.3), that is the order in approximating b and b% 

by fe( n ) and b% ; , respectively. We consider the case where T is the equilat- 
eral triangle II3 with vertices at 1, e 2l?r//3 , and e 4l7r / 3 . Then, by using the 
Schwarz-Christoffel formula it is not difficult to see that the coefficients b n 
of the associated conformal map (1.3) are given by bo = and 

cap^X-ir^ 3 )!, if 1 = 1, (41) 
0, if 1, { ' ' 

for n = 3m - l,m£N and / G {0, 1, 2}, where ( 2 ^ 3 ) denotes the binomial co- 
efficient; see, e.g., [13, §7.8]. Furthermore, it follows by using the properties 
of hyper geometric functions that 

3 rfl/3) 3 

b = cap(n 3 ) = -r^- = 0.730 499 243 103 •• • , (4.2) 

2 47H 

where T{x) denotes the Gamma function with argument x. 

By using the rotational property of the equilateral triangle, as this is 
reflected in the relation (2.11), it is easy to see that 

6 n _ M = 0, if k f {2,5,8,...}. 

This is actually the reason why we consider the two approximations b^ and 
br^\ Accordingly, we let and denote the two errors 

t (n). = b _ b (n) and t ("). = 62 _ 6 ("). (4.3) 

Then, from Theorem 2.1 we have that 

\t (n) \ < c{T)- and |4 n) |<c(r) * n G N. (4.4) 



n yn 
In Tables 4.1 and 4.2 we report the computed values of b^ n \ and 

(n) (n) 

&2 , ^ > w hh n varying from 100 to 200. We also report the values of the 
parameter s, which is designed to test the two hypotheses 

\t^\^l/n s and \t { ™ ] \ » l/n s . 

This was done by estimating s by means of the two formulae 

s n ■= log f ,,!* L I / log ( ^TT^l and 




| t (n+io)| 1 b ^ 10 



In view of Remark 3.1, regarding the stability properties of the Arnoldi 
GS process, we expect all the figures quoted in the tables to be correct. 
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n 




tin) 


s 


100 


0.730 487 539 


1 . 1 1 e- 


05 


1.9627 


110 


0.730 489 536 


9.70e- 


06 


1.9659 


120 


0.730 491062 


8.18e- 


06 


1.9685 


130 


0.730 492 255 


6.98e- 


06 


1.9708 


140 


0.730 493 204 


6.03e- 


06 


1.9728 


150 


0.730 493 973 


5.26e- 


06 


1.9745 


160 


0.730 494 603 


4.63e- 


06 


1.9761 


170 


0.730 495127 


4.11e- 


06 


1.9774 


180 


0.730 495 567 


3.67e- 


06 


1.9786 


190 


0.730 495 940 


3.30e- 


06 


1.9799 


200 


0.730 496 259 


2.98e- 


06 





Table 4.1. Equilateral triangle: Errors and rates in approx- 
imating b = 0.730499 243 103 •• • by M n ). 



n 


b in) 


An) 


s 


100 





243 555 903 


-5 


61e- 


05 


1 


9873 


110 





243 546 213 


-4 


64e- 


05 


1 


9886 


120 





243 538 830 


-3 


90e- 


05 


1 


9897 


130 





243 533 076 


-3 


33e- 


05 


1 


9907 


140 





243 528 504 


-2 


87e- 


05 


1 


9914 


150 





243 524 812 


-2 


50e- 


05 


1 


9921 


160 





243 521 788 


-2 


20e- 


05 


1 


9926 


170 





243 519 280 


-1 


95e- 


05 


1 


9931 


180 





243 517177 


-1 


74e- 


05 


1 


9936 


190 





243 515 396 


-1 


56e- 


05 


1 


9939 


200 





243 513 875 


-1 


41e- 


05 







Table 4.2. Equilateral triangle: Errors and rates in approx- 
imating b 2 := 0.243 499 747 701 •• • by . 



It is 



interesting to note the following regarding the presented results: 

The values of decay monotonically to b. 

The values of b^ increase monotonically b 2 - 

The values of the parameter s indicate clearly that 

\t^\ « 1/n 2 and \t^\ « 1/n 2 , 
This suggests that the two estimates 

\ t in)\ < c (r)I and |4 n) | < c(r)4=, n G N, 
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predicted by Theorem 2.1 are pessimistic. 

5. Proofs 

5.1. Proof of Theorem 2.1. The derivation in the case where T is piece- 
wise analytic without cusps is based on results from [12] and [13]. In par- 
ticular, we utilize the following fact (see [12, Thm 1.1]): 

n + 1 7 2(n+l) 

To — = 1 - On, (5.1) 

7T A n 

where 

0<a n < Cl (T)-. (5.2) 
n 

We also note the following estimate [12, Thm 1.2]: 



where 



j!l±lq,n {z) ^ {z){l + An{z)h (53) 

\A n (z)\ < - . ^L,, m 4= + c 3 (r) -. (5.4) 
1 Wl ~ dist(z,T) \&(z)\ yjn y 'n v ' 

Recall that we use c k (T), = 1,2,..., to denote a non-negative constant 
that depends on V only. 

The result (2.2) follows immediately from (1.9) and (5.1)-(5.2). For the 
general case (2.3), our proof relies on the use of the auxiliary polynomial 

yi+l 

q n -i(z) := G n {z) - — — p n (z), neH. (5.5) 
An 

This is a polynomial of degree at most n — 1, but it can be identically zero, 
as the special case when G is a disk shows. 

Next we fix k = 0, 1, 2, . . .. Then from (1.21), in conjunction with (5.5) 
and the orthogonality of p n , we deduce for any n = k, k + 1, k + 2, . . ., that 

— : K~k,n = {z— p n ,Pn-k) = (zG n - Zq n -i,p n _ k ) 

A« An 

= {zG n ,p n - k ) - (zq n -l,Pn-k) 

k 

= b(G n+ i,p n - k ) + ^2bj(G n -j,p n - k ) - (zq„-i,p n ~ k ). (5.6) 

3=0 

Thus, it remains to estimate the two different types of inner products 
appearing in (5.6), namely (pi,G m ) and (zq m ,pi). This is the objective of 
the following two lemmas. 

Lemma 5.1. Assume that T is piecewise analytic without cusps. Then, for 
I = 0, 1, 2, . . ., it holds that 

<"•<*■>={ r v *" :::!'+u + 2,..., <"> 
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where 

\U\ < ci(r)-. (5.8) 
m 

Proof. For the special case where m = I the result is a trivial consequence 
of the orthonormality property of the polynomial p m and the fact that G m 
is a polynomial of exact degree m with leading coefficient 7 m+1 . That is, 



{Pm, G m ) = ( Pm , 7 m+1 Z m + •••) = { Pm ,7 m+1 Z m ) 



l m+1 (Pm,T-Pn 



7 



m+l 



A, 



(5.9) 



Assume now that m £ {I + 1, 1 + 2, . . .}. Then, an application of Green's 
formula, the splitting (1.18) and the residue theorem give: 



(Pi,G r , 



Pi(z)G m {z)dA(z) 



G 



F ™ +l{z) 'dA(z) 



Pi\ z , 
G m + 1 



7T f 1 

m + 1 \ 27ri 7 r 

— (-/ 

m + l [ 27rz 



^(z)F m+ i(z)dz 



pi{z)<5> m+1 (z)dz + 



2ttz 



pi(z)E m+ i(z)dz 



7T 



1 



m + l ^2vri 7 r $ m+1 (z) 
1 



dz + ^~f Pl(z)E m+1 (z)dz^ 



pi(z)E m+ i(z)dz. 



(5.10) 



2(m + l)i 7 r J 

To conclude the proof we use the estimate given in [13, Lem. 2.5], to 



obtain 



[z)dz 



< Q2(r)||pi|| L a (G ) 



/ \E' m+l {z)\ 2 dA( 



z) 



1/2 



, (5.11) 



where we made use of the fact that E' m+1 £ L 2 (f}) (see Lemma 1.1) and 
that a piecewise analytic without cusps Jordan curve is quasiconformal and 
rectifiable. 

Therefore, from (5.10), the second relation in (1.20) and (5.11), we have 

1/2 

\H m {z)\ 2 dA{z) 



\{vi,G m )\ < C3 (r) 



(5.12) 



and the required result follows, because the last integral is 0(l/m 2 ); see [13, 
Thm2.4]. □ 

Lemma 5.2. Assume that V is piecewise analytic without cusps. Then, for 
every m £ N and I = 0, 1, 2, . . ., it holds that 



\(zq m ,pi)\ < ci(r) — . 

m 



(5.13) 
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Proof. The result is a simple consequence of Corollary 2.1 in [13] which 
states 

lkmllz,2(G) < c 2 (r)— , 

and the Cauchy-Schwarz inequality: 

\(zq m ,Pi)\ < lkgm||i2(G)lbdU 2 (G) < max{|z| : z G T} ||<? m || L 2 {G) . 

□ 



Returning to the proof of Theorem 2.1, we apply the results of the two 
previous lemmas to (5.6) and use (1.23) to obtain: 

n+l *pl 7 n-fe+l 

-£>n-fc,n = b(Gn+l,Pn-k) + 2_^bj(G n -j ,p n _ k ) + h~ (zQn-l,Pn-k) 



A, 



j=0 



fc-i 



' V ( ra - i) i 



1+0) 



7 



n-fc+l 



A 



n—k 



(5.14) 



where < w < 2, and O does not depend on n or fc. Furthermore, from 
(5.1)-(5.2) we have: 



.,n+l 



A„ 



7T 



n + l 



and 



A 



n—k 



n-k+l 



7' 



n — A; + 1 



7T 



1 + 



l + O 



J? 



n — k + 1 



(5.15) 



(5.16) 



Thus, by multiplying both sides of (5.14) by \ n ~k/l n k+1 we get 



K-k 7 



n+l 



^n—k,n &k 



A 



n—k 



7 



n—k+l 



V 7 .7=1 



(n - j) 



which, in view of the estimates (5.15)-(5.16), yields for n > k > 0, n > 1, 
that 



n — A: + 1 
n+l 



K-k,n = bk 



1 + 



1 



+0(Vn-fc + l) 



n- k + 1 

k-l 



\nj \{n- 3 ) 3 l +") 



(5.17) 



where an empty sum equals zero. This leads, for fixed k and n —> oo, to the 
required estimate (2.3), where now O depends on k but not on n. □ 
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5.2. Proof of Theorem 2.2. If T G C(p + a), with p + a > 1/2, then the 
following asymptotic formulas hold as n — > oo, see [14, p. 19-20]: 

n + 1 7 n+1 _ / 1 \ , 

i + obra ( 5 - 18 ) 



7r A n \^n 2 (P+ a ) 
and 



a/^^(^)(i + o(^) 



7T 



The proof of the theorem goes along similar lines as the proof of Theo- 
rem 2.1 given above. More precisely, for deriving the result for 6 n +i in we 
use the estimate (5.18) in the place of (5.1)-(5.2). 

For the general case k = 0, 1, . . ., we need estimates for the inner products 
{pi, G m ) and {zq m ,pi}. This is done in the following two lemmas, which play 
the role of Lemma 5.1 and Lemma 5.2 in the proof of Theorem 2.1. 

Lemma 5.3. Assume that T G C(p + l,a), with p + a > 1/2, then for 
I = 0, 1, 2, . . ., it holds that 

<"•"->-{ eT ,/A "' ™ = *+m + 2,..„ ( 5 - 2 °> 

where 

^^c^T)-^^. (5.21) 

Proof. The result for m = I is established in Lemma 5.1. Hence, we only 

consider the case m = I + 1, 1 + 2, 

The following estimate has been obtained by Suetin for T G C(p + a); see 
[14, Lem. 1.5]: 



— / H m {z)E m+ i{z)dz 
2m J r 



^ 2 ( r )^y- (5-22) 

By using Green's formula in the unbounded domain 0, together with (1.20), 
it is readily seen that 

/ H m {z)E m+l (z)dz = _^±1 / IfZ^JpdAO*). ( 5 - 23 ) 
Hence, from (5.22), 

J a \Hm{z)?dA{z) < C3 (r) m2(p ' Q)+1 ) (5.24) 

and the result (5.21) follows from the estimate (5.12), which is applicable 
in this case because any smooth Jordan curve is also quasiconformal and 
rectifiable. □ 

Lemma 5.4. Assume that T G C(p+ l,a). Then for every m G N and 
I = 0, 1, 2, . . it holds that 
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Proof. As in the proof of Lemma 5.2 we have 

\(zq m ,Pi}\ < max{|z| : z G T} \\q m \\ L 2( G y 
The result of the lemma then follows from (5.24) and the estimate 

1/2 

\q m \\L2(G) < c 2(r) / \H m {z)\ 2 dA(z) 

established in [13, Thm 2.1] for domains bounded by a quasiconformal and 
rectifiable boundary. □ 

In order to conclude the proof of the theorem, we need an estimate for 
the decay of the coefficients b n , when the boundary T belongs to the class 
C{p + l,a), with p + a > 1/2. This is done in [13, Cor. 1.1], where it is 
shown that 

IM<c 3 (r)-^l^, iiGN. (5.26) 

Therefore, by using the results for (pi,G m ) and (zq m ,pi), obtained in the 
previous two lemmas, together with (5.18) and (5.6), we see that 

fe-l 



1 



n+l 



A, 



'n—k,n 



o 



1 



n p+a+l/2 



1 



{j{n _ j))P +a+l/2 



1 



n-k+l 



A 



n—k 



where O does not depend on n or k. Furthermore, from (5.18) we get 

,,n+l 



7 



A, 



7T 



n + l 



1 + 



1 



n 2(p+a) 



and 



A 



n—k 



n-k+l 



1 



n — k + 1 



7T 



1 + 



(n-k + 1) 2 (p+ q ) 



(5.27) 



(5.5 



The above yield, for n > k > 0, n > l, that 

1 + 



n — k + 1 



n + l 



n—k,n 



1 



o 



( n p+a+l/2 ) + 

j=l 



(n-k + 1) 2 (P+ Q ) 
fe-i 



1 



(i(n _ i))P+ a+l/2 



+ 0(^-/^ + 1) 

(5.29) 



where a empty sum equals zero. This leads, for fixed k and n — > oo, to the 
required estimate (2.5), where now O depends on k but not on n. □ 

5.3. Proof of Theorem 2.3. Assume that T := dG is an analytic Jordan 
curve. Then the conformal map $ has an analytic and univalent continuation 
across r in G. Let g < 1 be defined by 

g := inf{r : $ is analytic and univalent in ext(L g ) \ oo}. 
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Then the following asymptotic formulas of Carleman [3] hold asn-s-oo: 



n + 1 7 n+1 

7T A,„ 



i + o{ef 



(5.30) 



and 



n + 1 



7T 



Pn{z) = 

see [6, p. 12]. In particular 
where 



$ n (z)&(z){l + 0(V^Q n )} , zen, (5.31) 



7 



(5.32) 



= ^ i/^w"- 1 -" - ^ i/a^-^- 1 "™, w = $(z), (5.33) 



with 



-£v\A v f + J2v\a„fg-^< 



J2n+2 



(n+l)(l- g 2n )' 



n G N; (5.34) 



see [6, p. 15]. 

We fix two different points z and zq on V and define 

Oi+iW := £pi(CK+ (j . + ^ +1 ^ +1 (^0). 

Then, by using integration by parts and the change of variable w = 3>(C), 
we have from (5.32) and (5.33) that, for any j G N, 



-0 



(j + 1)7 J+1 
(J + 1)7 J+1 



& +1 {z) + 



A, 



■j'+i 



ii)0 



,l/=l 



I/=l 



where u>o = 3>(zo)- We claim that for = 1 there holds 

3 oo 



v=l v=l 



W 



o 



' logp' + i) 



(5.35) 
(5.36) 



Indeed, 

3 



^2 A v w v + ^ a v w~ v 

u=l v=\ 



i/=i i/=i 
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< 



3 1 



\ u=l \ V=\ \ V=\ 



v\a v \ Q 



2 n -2u 



g zu 



< Cl (r)Viog(i + i) 



^>|A,| 2 + c 2 (r), 



< C3(r)i 



i + i 



(5.37) 



by (5.34), which establishes the claim. 
Hence, using the estimate (5.30) we get 

Qi+i CO " 



= {i + o(V(i + i)iog(j + i))^'}, ^er. (5.38) 



Next, by Green's formula we have for fixed k = 0, 1, . . . and re > fe + 1: 



2vri 



n — + 1 . 
n + 1 



-*n— fc,n 



n — k + 1 

2m \ r~i \ZPn,Pn-k) 

V n + 1 
27ri In — k + 1 
~2i \ n+l 

j z® n (z)&(z)$ n - k+1 (z)dz + /t, 

$ n (z) 



zp n (z)Q n ^ k+1 (z)dz 



r $"- fc +!(z 



u,l=l w 

27ri6 fc + /in, 



<&'(z)zdz + /i n 



(5.39) 



where 



/i n = O^f?" + O (v> - k + 1) log(n - k + 1)J + 0(V^)Q n }. 

(5.40) 

Thus, for > fixed and g < 1, 



n — A; + 1 , 
n + l 



-*n— fc.n 



6^ + 0(y / n logng n ), as n — >• oo. 



(5.41) 



It remains to prove (2.6). This follows at once from the strong asymptotics 
for the leading coefficient (5.30) and the relation (1.9). □ 

5.4. Proof of Theorem 2.4. We first note that our assumption (2.8), com- 
bined with (1.9), implies that 



lim sup 

n— >-oo 



n + 2 A n 
n + l An+i 



l/n 



< 1. 



(5.42) 
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Now set 

£ •= j n + 2 Xn 1 1 

V n + 1 A n+ i 6 

so that 

limsup|£n| 1/n < 1. (5.43) 

n— >oo 

At the other hand, we have from (5.1) and (1.4) that 

(1 + £n) 



1 - "n+1 



1 — a. 
Hence, 

a n - a n+ i 



s " (l-a n )(2 + £ n )' 

and by using the fact that £ n — > 0, as n — > oo together with < a n < 1 and 
a n — > 0, as n — > oo, we obtain the double inequality 

ci|a n - a n+ i\ < |£ n | < c 2 \a n - a„+i|, (5.44) 

for some positive constants c\ and c 2 . 

Now, by expanding a n in the telescoping series 

a n = («n - «n+l) + («n+l ~ «n+2) H , 

we conclude, in view of (5.43)-(5.44), that 

limsupay™ < 1, (5.45) 

n— >oo 

and this, in view of Theorem 1.3 in [13] leads to 

limsup l^nl 1 /™ < 1. 

n— >oo 

The last inequality implies that the conformal map *$>(w) has an analytic 
continuation across T into B (see (1-3)) and thus T is the analytic image 
of T. Therefore, around any wq G T, the map ^ can be represented by a 
Taylor series expansion of the form 

*S>(w) = ^(w ) + ai(w - w ) + a 2 (w - wq) 2 + a 3 (z - z ) 3 

If we had $>'(w ) = 0, then 

#(u>) = *(wo) + a 2 (w - w ) 2 H , 

with a 2 7^ 0, because ^ is univalent in A. These show that wq would be 
mapped by ^ onto an exterior pointing cusp on T. Since, by assumption, 
this cannot happen, we see that ^?'(w) ^ 0, w £ T, which yields the required 
property that T is an analytic Jordan curve. □ 
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5.5. Proof of Theorem 3.1. Recall that n := 2m. On \w\ = R, where 
q < 1 < R < oo, we have from (1.3) and (3.2) 

m ir(n) ,i oo I , I 

i^h-^ ) wi<i^ ) -^+E^f M + E l §- 

k=0 k=m+l 

Therefore, by using the result of Theorem 2.3 (see also (5.39)) and the 
estimate 

k 

\b k \ < Cl (T)^=, kN; 

see [13, Cor. 1.1] we get 

|^H-^HI < C2 (r)/ m i? + C3 (r)y^no^^ + C 4(r)(|) m , (5.46) 

which yields the desired estimate. □ 
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